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Abstract. We present the methodology and some preliminary results of our study 
of the relationship between a Mira's pulsating photosphere and its surrounding molec- 
ular layer(s) throughout several pulsation cycles, based on spatially resolved data. Our 
dataset consists of archival narrow-band observations in the near-infrared H and K 
bands obtained with the Palomar Testbed Interferometer between 1999 and 2006, ex- 
tended with a few nights of VLTI AMBER low spectral resolution data and near- 
infrared SAAO photometry. The fitted model is the geometric star + layer model pro- 
posed by Perrin et al. (2004), in which the physical parameters (diameter and temper- 
ature of star and layer; wavelength dependent optical depth of the layer) are given a 
sinusoidal time dependence. 



1. Introduction 

The issue of Mira pu lsation modes has l ong been an issue of diameters. Ve locity ampli - 



tude measurements dBessell et al.lll996l) and period-luminosity relations ( Woodll2000h 



are eviden ce in favor of the f undamental pulsation mode hypothesis, while measured 
diameters dTuthill et al.1 11994) in the past pointed towards a first overtone pulsation. 
Moreover, the scales of the varying Mira diameters with wavelength can not be ac- 
counted for by the trend in continuum opac ity reaching a minimum arou nd 1.6 fxm. 
This changed when iMennesson et al . (2002) (and later iPerrin et al.l 120041) found that 
the measured diameters are actually biased by the presence of a molecular layer (con- 
sisting mainly of CO and/or H2O) at up to a few stellar radii above the photosphere. 
The observed strong diameter variations with wavelength are in this simple geometri- 
cal model explained by the varying transparency of t he sh ell and the wavelength de- 



pendence of black body radiation. iLe Bouquin et all d2009h later confirmed this basic 



picture by performing a model-independent image reconstruction on an interferometric 
dataset continuous in both the spectral (1.5 - 2.4 /urn) and spatial domain for the Mira 
variable T Lep. 

Although the general picture is now well established, a lot of questions related 
to the molecular layer remain unresol ved. What is its role in driving the stellar wind? 
What are its exact sources of opacity? Le Bouqui n et all ([2009) needed a source of con- 



tinuum opacity (r = 0.2) on top of CO and H2O to explain the optical depths required 
by the fit to the interferometric observables. How do the photosphere and molecular 
layer vary with the pulsation? How does this time dependence compare with predic- 
tions of pulsation theory? Are Miras really pulsating in the fundamental mode? 
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The immediate goal of our study is to quantify the temporal behavior of photo- 
sphere and molecul ar layer for a samp le of Mira variables. To this end, we extend the 
analytical model of IPerrin et all (120041) by assuming a sinusoidal time dependence for 
the physical parameters 9+, 7\, 9i ayer , Ti ayer and t(A). 



2. Data and finding the best fit 

Due to its large and diverse sample (~ 100) of observed Mira variables, and excellent 
temporal coverage (about 8 yrs), we tried to apply our model to the publicly released 
Palomar Testbed Interferometer (PTI) dataset. The PTI is a fixed three telescope in- 
terferometer (baselines of 86m, 87m, and 110m roughly oriented NW, SW, and NS 
respectively) operating in five narrow spectral bands covering 2.0 - 2.4 nm or (i.e. not 
within the same night) four channels covering 1.50 - 1.75 fitn. 

R Vir is the second best observed Mira in the PTI archive and has about 275 



squared visibilities per wavelength, spread over 12 pulsation cycles. Eisner et "ai1 ( 20071) 



previously applied the same model to a spectrally dispersed interferometric Keck dataset 
of R Vir during rise to maximum visual light. The Keck data, which have a higher spec- 
tral resolution (R ~ 230 vs. R ~ 30) but carry less spatial information, are consistent 
with the PTI data at similar phases. We extended our dataset with three VLTI-AMBER 
points in low spectral resolution on triplet A0-G1-K0 during GTO time in February 
2009. In contrast to what was presented on the conference, there is no discrepancy 
between these AMBER data and PTI data at similar phases (the presented discrepancy 
was caused by a local software problem in the calibration process). 

To fit the data we make use of two independent methods: a standard Nelder-Mead 
Simplex algorithm and a genetic algorithm. Our main concern at the moment are the de- 
generacies the 2-component model gives: very similar visibility curves can be produced 
with very different values for the physical parameters. From our fitting and simulations 
we conclude that the PTI dataset has insufficient spatial information at high enough 
accuracy to constrain the model. The AMBER data do constrain the model, but only at 
that particular phase. Hopefully, combining the two datasets gives enough constraints 
to disentangle the temporal behavior of the different components in the model. 
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